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We study a two-dimensional electron system in a magnetic 
field with a fermion hardcore interaction and without disor- 
der. Projecting the Hamiltonian onto the n th Landau level, 
we show that the Hartree-Fock theory is exact in the limit 
n — » oo, for the high temperature, uniform density phase of 
an infinite system; for a finite-size system, it is exact at all 
temperatures. In addition, we show that a charge-density 
wave arises below a transition temperature Tt. Using Landau 
theory, we construct a phase diagram which contains both 
unidirectional and triangular charge-density wave phases. We 
discuss the unidirectional charge- density wave at zero tem- 
perature and argue that quantum fluctuations are unimpor- 
tant in the large n limit. Finally, we discuss the accuracy of 
the Hartree-Fock approximation for potentials with a nonzero 
range such as the Coulomb interaction. 



I. INTRODUCTION 

In a two-dimensional electron system (2DES) in a 
strong magnetic field, the fractional quantum Hall ef- 
fect [ |IJ may be observed when electrons partially fill the 
lowest (n = 0) Landau level. In this paper, we consider 
what happens if it is a high Landau level that is par- 
tially occupied. The large-n limit generates a problem 
that we can, in part, solve exactly. To this extent, the 
high Landau level limit is analogous to the limit of high 
dimensions for spin models [ || or the Hubbard model [ 

1- 

Electron states in lower excited Landau levels (n = 
1,2) have received some attention in the past [ ||], and 
recently Aleiner and Glazman [ |^] have derived an ef- 
fective interaction for electrons in high Landau levels, 
taking into account the wavevector dependent dielectric 
constant due to screening by lower, fully occupied levels. 
Using this interaction, Foglcr, Koulakov and Shklovskii 
(FKS) [ D have studied a 2DES in a weak magnetic 
field at zero temperature, within the framework of the 
Hartree-Fock approximation. For a filling fraction, v, of 
the partially filled Landau level near v — 1/2, they find 
a phase of fully occupied (incompressible) regions alter- 
nating with regions of unoccupied orbitals. These regions 
are translationally invariant in one dimension ('stripes'). 
At lower filling fractions, these stripes break up and there 
is a phase of fully occupied islands ('bubbles') on a tri- 
angular lattice in an otherwise empty system. The width 
of the stripes and the lattice constant are each of order 
of the Larmor radius. Experimental consequences of this 



structure, in particular for tunneling into the 2DES, have 
been investigated extensively by FKS. 

Both the interest and the difficulty of studying elec- 
tron correlations in a partially filled Landau level stem 
from the fact that the single-particle problem is macro- 
scopically degenerate. Interactions set the only impor- 
tant energy scale and the problem is intrinsically non- 
perturbative. In the n Landau level, however, a small 
parameter does appear, which manifests itself in the 
existence of t wo new lengthscales besides the magnetic 
length, I = y/hj (eB): Af = 2-K/kp = 2nl/y/n, the Fermi 
wavelength, and the Larmor radius, Rl ~ l\/2n, the size 
of the smallest wavepacket that can be constructed in the 
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Landau level. We use this small parameter, 



, to 

obtain a problem that is amenable to exact solution. 

We consider a 2DES assuming that the partially filled 
Landau level is spin-polarized Using a fermion hard- 
core potential of strength v -C huj c as a model interac- 
tion between the electrons, we are able to resum per- 
turbation theory at high temperatures, for n — > oo. 
(hu> c — heB/m is the cyclotron energy and from hereon, 
we set I = h = 1). 

We find that, for the high-temperature phase, the 
Hartree -Fock approximation is exact in the large-n limit 
(section III B ); for a system on a cylinder of finite circum- 
fere nce, it remains exact down to zero temperature (sec- 
tion [II A ) and the eigenstates are single Slater determi- 



nants of Landau orbitals. The high temperature phase is 
unstable towards the formation of a charge-density wave 
(CDW) of wavelength A min ~ 2.9R L (section [HIPP - We 
construct a phase diagram for the system using a Landau 
expansion (section IV). We find that there exist both 
a unidirectional CDW, near v = 1/2, and a triangular 
CDW away from v = 1/2. 

Our work hence generalises to finite temperatures the 
results obtained by FKS at zero temperature. It also 
represents a generalisation of the study of charge-density 
waves in the lowest Landau level by Fukuyama, Platz- 
mann and Anderson [|tJ. In particular, we show that the 
transition from the uniform state near v = 1/2 is to a 
unidirectional rather than triangular CDW. 

We also (section VI) compare the theories with hard- 
core and finite range potentials. It emerges that the 
Hartree-Fock approximation is not exact for potentials 



of finite range R (section VI A ) even in the large- n limit 



We show, however, that corrections to the Hartree-Fock 
approximation are governed by the parameter R/l, which 
is small in the weak-field limit investigated by FKS. For 
potentials where R/l is not small, we demonstrate ex- 
plicitly that there exist uniform density states which, for 
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these interactions, have a cohesive energy of the same 
order as that of the CDW states. It therefore remains 
possible that the groundstate in high Landau levels may 
depend on details of the interaction, and may not neces- 
sarily have CDW order. 



II. THE HAMILTONIAN 

We start by considering the Hamiltonian projected 
onto one spin orientation for the n th Landau level: 

H = - ^2 {rq\V\pk)ala\a k a p 

r,q,p,k 



k,m ' 



6 (<ly ~ l)a{ +l ai +m a k+m+l a k 

M(m,l) al +l al +m a k+m+l a k . 



(2.1) 
(2.2) 



Here, creates an electron in the Landau orbital (r\p) — 
tpp (r) of the n th Landau level with pseudomomentum p, 

tpp(r) = L~^e wx cj)^(x — p), where is the n th har- 
monic oscillator wavefunction. Periodic boundary condi- 
tions are applied in the y-direction for a system size L. 
Also, a n (q 2 /2) = [u n {q 2 /2)) 2 = (L n (q 2 /2) e -9 2 /4)2 ; Ln 
being the n th La guerre polynomial. The Fourier trans- 
form of the potential is V"(q) = / cPr e lq r V(r), and 
all quantities / with a tilde are to be understood as the 
Fourier transform of the function /. To derive (2.1), we 
have used f HI 



(fc'|e"i- r |fc) = ^S(k' 



(2.3) 



We represent the vertex M(m, I) , which is independent 
of k, as in Fig. [|. 

In this paper, we consider primarily the simplest pos- 
sible potential: a contact interaction for fermions. In the 
lowest Landau level this interaction is known to repro- 
duce much of the physics of the fractional quantum Hall 
effect [ ^| and the Laughlin wavefunction [ [l(J is known 
to be its exact [ [ll| ground state at v — 1/3. It is highly 
singular at the origin in real space but has a well-behaved 
Fourier transform 



Vkc(g) = -vq 2 



(2.4) 



v being the interaction strength. 

In addition, in section VI we discuss the screened 
Coulomb interaction, Vagi recently derived by Aleiner 
and Glazman [ M , as well as two simple finite range po- 



tentials, V\ and V2. The potential Vag (2.5) is an effec- 
tive intra-Landau level interaction potential which takes 



into account the presence of lower, fully-occupied Lan- 
dau levels through the dielectric constant cag- It has the 
following important features: there is little screening for 
q ^5 n^ 1 and q > 1, but for q ~ n~5, eag ~ Rl/gb- 
In weak magnetic fields, this is very large, of order n', 
whereas in strong fields, 6ag ~ 1 for all q. The explicit 
form of Vag is 



V AG (q) = 



2Tre 2 



47re Ke A G (?)<?' 



where 



£Aa(q) = 1 



qa B 



[1 - J 2 ( 9 R L 



(2.5) 



(2.6) 



Jo is the Bessel function, 0,3 = ti 2 Aneon/ (me 2 ) is the 
Bohr radius, and k is the relative permittivity. 

The two finite-range potentials we consider are of range 
R. One potential is box-like, non-zero and constant up 
to a distance R and zero for larger distances; the other 
one decays exponentially with a decay length R. We 
take Af <C R <C Rl- The Fourier transforms of these 
potentials are 



Vi(p) =v 1 R 2 J 1 (qR)/qR- 



V 2 (p)=v 2 R 2 (l + q 2 R 2 )-i 
(2.7) 



where the Vi are the strengths. Potentials of a finite range 
arise, for example, if there is a metallic gate close to the 
two-dimensional electron system which effectively screens 
the interaction at distances larger than the 2DES-gate 
separation. 



III. THE LIMIT n -> 00 

A. The vertex M(m,l) 

We first investigate qualitatively the behaviour of 
the individual matrix elements. To do so, we replace 
a n (q 2 /2) by its WKB envelope [ [l^], normalized to 
be consistent with J Q (Ltq(x)) 2 e~ x dx = 1. We set 
a n (q 2 /2) = for q 2 > 8n, and for q 2 < 8n 



My) 



1 



7T ^/q 2 (8n - q 2 ) 



(3.1) 



This is an asymptotic expression away from q = and 
q = \/8n. In the vicinity of these points, we have used 
the following expressions [ [ll| : 



My) = JoiQ^n) 



Mv) = (Ho) 183 ) 



q 2 < n (3.2) 
! -8n<n* (3.3) 
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In considering the envelope only, we neglect the rapid 
oscillations of a n (q 2 /2), since the potential does not have 
structure on this scale. 

To estimate M(m, I) for the hardcore potential we 
rescale q x = \/8n Q x and find 



M(m,l) 



^2 
L jo 



^V^l(-v(Ql + f )) x 



cos(Q x mv8n) 



(3.4) 



For large n, the oscillating cosine suppresses all 
M(to, I) for which to 3> n~ 3 . 

The consequences of these simplifications are particu- 
larly clear if we consider a system on a cylinder of finite 
circumference L. Then, m = 2?:a/L, with a integer, so 
that, in the limit n — ► oo with L fixed, only M(m = 0, 1) 
is non-zero [ [l4| . The Hamiltonian in this limit thus con- 
sists only of exchange terms (to = 0). It is therefore diag- 
onal in a basis of Slater determinants of Landau orbitals. 
The charge density of any single such Slater determinant 
is independent of the y-coordinate. The problem of find- 
ing the ground state is therefore reduced to finding the 
unidirectional charge-density wave which is lowest in en- 
ergy, and the Hartree-Fock approximation for this system 
is exact. 

Note that the spatial dependence of the exchange in- 
teraction is as follows: M(m = 0, 1) varies little as a 
function of I for |Z| < 2Rl- For |/| > 2Rl, the cyclotron 
orbits do not overlap. The effective range of the contact 
interaction is therefore determined by the spatial extent 
of the wavefunctions and is thus of order of the Larmor 
radius. 

It is interesting to compare our conclusions with those 
of Rezayi and Haldane [ who considered a similar 
problem in the lowest Landau level. They found by nu- 
merical diagonalisation that a single Slater determinant 
of Landau orbitals is the ground state if the magnetic 
length is much larger than the circumference of the cylin- 
der. In the system we consider, the important lengthscale 
is Rl and the limiting behaviour is reached for L <C Rl- 



B. Diagrammatic analysis 

We now turn to the infinite system and show by means 
of a diagrammatic analysis that, in the absence of sym- 
metry breaking, the Hartree-Fock approximation is exact 
for the hardcore potential, in the high temperature phase, 
in the limit n — > oo. To do this, we group all diagrams 
for the self-energy into three classes (Fig. g): first, those 
diagrams which contain neither closed fermion loops nor 
crossed interaction lines (class A); second those which do 
not contain closed fermion loops but do contain crossed 
interaction lines (class B); and third those which contain 
closed fermion loops (class C). 



For the analysis of the n-dependence of these diagrams, 
we note that the sum over Matsubara frequencies and the 
integral over intermediate states decouple: the bare prop- 
agators, k , which are diagonal in the pseudomomen- 
tum k, are the same for all states in a given Landau level, 
since G^ k = (iuj m ~ (e n (k) - /i))" 1 = (iu) m - (e-^))" 1 . 
There is no n-dependence in the propagators since the 
dependence on n in e„ is absorbed into an n-dependent 
chemical potential /i. Therefore, only the behaviour of 
the integrals over intermediate states has to be consid- 
ered when taking the large- n limit. 

For a diagram in the first class with one interaction 
line, the integral over intermediate states is 

^ Svn f 1 Q 2 dQ 2vn 

> M(0,q y ) = T l = 3.5 

~ 1" JO \/l - Q 2 7T 



This class of diagrams in fact makes the leading con- 
tribution to the self energy for n — > oo and hence we 
scale v with n so that vn is constant. Next wc sum all 
the diagrams in this class; in appendix |X| we establish an 
upper bound of (vn) x n~3 to the contribution from other 
diagrams, which vanishes for n — > oo. The sole exception 
are the diagrams which contain only Hartree loops (Fig. 
||) and interaction lines which do not cross; these van- 
ish for our potential since Vhc(0) = 0. If this were not 
the case, we would have to compare their contribution 
with that of the Fock diagrams when considering what 
to retain for large n. 



C. The self energy Efc 

We sum the diagrams which do not vanish as n — > 
oo using the Dyson equation (Fig. ||), noting that the 
absence of a Hartree contribution reduces it to the one 
in Fig. |5| so that 



l,m 



-M(0,l) 



„i-E„,0 + 



E 



/3 



E 



iE 



ni 



(e-M) 



p \_m,n r 



-M (0, m) 



-1 



iE~. 



l- /Mi- ^(0,0 



(3.6) 



For the hardcore potential, £ i= _ 00 M(0, 1) is negative 
(Eq. 3.5), and the self-energy has no poles. 

We see that is independent of k, because of the 
translational invariance of the system; to keep the filling 
fraction v fixed, Sfc will be compensated by an equal shift 
in the chemical potential. Hence, to leading order in n, 
the one-particle propagator remains unchanged. 
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D. The Bethe-Salpeter equation 

The two-particle vertex T, on the other hand, is af- 
fected by the interaction. An analysis like the one above 
applies also to the perturbation expansion of T, and in 
the end we only have to sum ladder diagrams illustrated 
in Fig. | [ 0- Tne Bethe-Salpeter equation (Fig. @), in 
the limit n — > oo, is 



i(uj + A) - (e - /i) i(A) — (e — //) 



T(q-l,u). 



Note that, because the constituent vertex M(m,l) and 
the propagators are independent of the pseudomomen- 
tum k and the Matsubara frequencies £1 and 5, so is F. 

For uj ^ the sum over A yields and T(q,uj) = 
M(Q,q), For uj = 

r(g,0)=M(0,«) + 

)9i/(l-i/)5^M(O,0r(g-Z,O) 



using second order perturbation theory with a Hartree- 
Fock Hamiltonian. These authors also show that, within 
Hartree Fock theory, a second order transition at T 2 is 
preempted by a first-order transition at a higher temper- 
ature T\ to a triangular CDW. 

In the following, we extend the work by FPA to the 
Landau level [ fL9[ . In the spirit of Landau theory, 
FPA expand the free energy in powers of the order pa- 
rameter S\ for the triangular phase using a diagrammatic 
calculation described in Refs. [ |2(]] and obtain 



SFx = NU{Q){a\S 1 \ 2 + b^ 3 + c^ 4 ) 
where the coefficients are given by 
a = 3(1 — v{\ — v)u) 
b = -2u 2 v{\ -v)\l- 2j/|cos(^-Q 2 ) 



(4.1) 



u 3 is(l - V) 



cos( 



6(i/ 



(4.2) 



and 



r( s ) 



M(0,s) 



1 - M(0,s)j3v(l - v) 



(3.7) 



where f and M are Fourier transforms of T(q, 0) and 
M(0, q) wit h respect to q. 2wM is in fact the Fock po- 
tential (5.5). 

The two particle vertex is divergent at a tempera- 
ture T = M{Q, s)v{\ — v). This indicates an instabil- 
ity towards the formation of a charge-density wave at a 
wavevector p m - la where M(0,p m i n ) is largest. We discuss 
the phase diagram of the charge-density waves in section 



[V . Had we chosen the potential to be attractive, the 
system would phase-separate into a fully occupied and 
an empty region below the critical temperature. 



IV. CHARGE-DENSITY WAVE INSTABILITY 
AND THE PHASE DIAGRAM 

In the preceeding section, we have shown that the uni- 
form system is unstable to formation of a CDW with a 
wavevector p mm at a temperature T 2 = v (l — v)U{p m \ xl ), 
with U(p) — uf(p)/(2tt). In section |VB| , we find for the 
hardcore interaction that p m in — 2.16/Rl- For a gen eral 
potential, replacing up by — uhf = -% + uf (5.4) in 
these and subsequent expressions yields a general analy- 
sis of the CDW phase diagram for a Hamiltonian in the 
Hartree-Fock approximation. 

The problem of charge-density waves in the lowest 
Landau level has been studied extensively J_ IS | , and 
Fukuyama, Platzmann and Anderson (FPA) [ [7| have ob- 
tained the same expression for the transition temperature 



Here, Q is the wavevector of the triangular phase de- 
termined below, u = U(Q)/T and N is the number of 
states in the system. By the same method, we find the 
free energy for the unidirectional CDW: 

6F 2 =NU(Q)(A\5 2 \ 2 + C\6 2 \ 4 ) 



where the coefficients are 
A = a/3 and C = u 3 v(l - v) 



1 1/ ^2 

- + -{v — Y 

8 2 V 2' 



(4.3) 



The temperature at which the first order transition 
occurs from the uniform phase to a triangular CDW is 



Ti =T 2 



U(Q) L , ecos 2 (^Q 2 ) 

U(p min ) y / + gcos 2(\/3Q2) 



(4.4) 



where e = - ±) 2 , / = i + §(„ - i) 



/' — i) 2 . Ti is found by maximizing (4.4) with respect 



3 v" 2 

to Q. For high Landau levels, this calculation simplifies 



as follows. Since Q 



Q 2 ) 



Hence, for large n, T\ is obtained by setting Q — p min . 
Therefore, T\ = T 2 (1 + e/(/ + g)) is always larger than 
T 2 , except for v — 1/2, where there is a triple point: T\ = 
T 2 = T t . In the remainder of the calculation, we can set 
cos 2 (-^Q 2 ) = 1 and U(Q) — U(p m in), thereby removing 
the second term in brackets in ( |4.2| ). 

We now turn to the question of which phase is present 
for T < T t . Since our expression for the free energy is 
valid only for small 5, we consider the neighbourhood of 
the triple point. Minimizing 5F\^ 2 with respect to 5\_ 2 , 
we find that the unidirectional CDW is more favourable 
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at v — 1/2 for T just below T t . As the filling frac- 
tion is lowered at constant T, there is a first-order phase 
transition to the triangular CDW. When v is lowered 
even further, the uniform density state becomes more 
favourable. We have also examined competition between 
the striped phase and a square CDW, using the expres- 
sion derived by Gerhardts and Kuramoto [ |l8); we find 
that at v = 1/2 for T just below T t , the square CDW 
has a higher free energy than the unidirectional one for 
n — > oo. The resulting phase diagram is displayed in 
Fig. H for v < 1/2. The results for v > 1/2 follow from 
particle-hole symmetry. The emergence of the unidirec- 
tional CDW can be explained from the expansion of the 
free energy as follows. As v approaches 1/2, b tends to 
zero. The gain of SFi relative to 8F2 due to its larger 
second-order coefficient, a = 3A, which is almost zero 
near T t , is offset by the cost incurred due to its larger 
fourth-order coefficient c > C, and the unidirectional 
CDW is energetically more favourable. 

As the temperature is lowered far below Ti, the or- 
der parameter grows and our expansion ceases to be 
valid. An additional difficulty is that higher harmon- 
ics will start to play a role: the density profile must 
be constructed by occupying states in the n th Landau 
level, and must in particular not require a filling fraction 
larger than one at any point. In section VB we show 



that the cost for higher harmonics is small. We therefore 
expect that these considerations alter only the shape of 
the CDW without changing its structure or periodicity. 
In this way, we expect competing unidirectional and tri- 
angular patterns to persist at wavevector p m i n down to 
zero temperature. 

If we extrapolate the phase-separation line T\% beyond 
the range of validity of our expansion to T = 0, we obtain 
a phase diagram (Fig. ^ with a first-order transition at 
v = 0.35 for T = 0. This is close to the value v — 0.39 ob- 
tained numerically by FKS, separating the unidirectional 
pattern from the triangular one. 



V. ZERO TEMPERATURE ANALYSIS 

A. The effective Hartree potential 

We start this section with a general remark concern- 
ing the Hartree-Fock approximation for a Hamiltonian 
projected onto a single Landau level. MacDonald and 
Girvin [ [2l]] have shown, using analyticity in the lowest 
Landau level, that the Fock potential can be expressed 
solely in terms of electron density as a function of posi- 
tion. This result is in fact not a consequence of the special 
analyticity properties of the wavefunctions in the lowest 
Landau level: if the Hamiltonian is projected onto any 
single Landau level, the Hartree Fock energy is simply 
1/2 J d 2 qV^.ff (q)|rl(q)| 2 , where n(q) is the Fourier trans- 
form of the electron density and the effective potential 
is 



^eff(q) = 

K F ff (q) = 



gmKq) 
1 



M<7 2 /2) 



/ 



ipq «4y) 



(5.1) 



(5.2) 



The crucial point is that after projection, the density 
matrix in the Landau orbital basis, p^^i can be e xpre ssed 
in terms of n(q). Details of the derivation of (5.2) are 
given in appendix 0. 



B. Self-consistent equations 



The diagrammatic analysis of section III B applies to 
the high temperature phase where the electron liquid has 
a uniform density. For that case, the Hartree-Fock ap- 
proximation is exact as n — > 00. We now allow for sym- 
metry breaking and solve the self-consistent set of equa- 
tions depicted in Fig. |]. For a system on a cylinder 
of finite circumference, we have shown that the Hartree- 
Fock eigenstates have a probability density that is trans- 
lationally invariant in the y-direction and hence, that it 
is sufficient to allow the self-energy to depend on one co- 
ordinate only. (The infinite system is discussed further 
in section |V C|). The Dyson equation reads 



2tt 



M(0,i-ft)- 



-M(Z-fc,0) 



E - Efc - (e - fi) 



B-E fc -(e-/i). 

= J2[M(l - fc,0) - M(0,l- k)](m) 
1 

where (rife) = — 0(£fc + (e — A*)) is the occupancy of the 
Landau orbital \k) and 8 is the Heaviside step function. 
This gives 



(5.3) 



where a n [p) \ ^'io^y) = -{vn) P 2 - , (5.4) 



-«f(p) = - / qdqV (q)a n ( — )J Q {pq) (5.5) 
^A(vn) 1 F 2 (^;l,2;-2np 2 ), (5.6) 

if 2 being a hypergeometric fu nctio n. 

N ote that the Fock potential ( |5.5| ) follows directly from 
( |5.2|) as it must. The factor of a n (p 2 /2) arises from the 
matrix elements of t he p otential between states in the n th 
Landau level. Eqs. (5.4) and ( fT5| ) were first obtained by 
FKS [ |6| who considered the interaction between guiding 
centres of coherent states in the n th Landau level. 

The asymptotic expression ( [5.6| ) for uf(p) holds as long 
as p < ^/n, and depends on p only via the product 
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pRl indicating that the important lengthscale for the ex- 
change potential is Rl- The functions uu(p) and —uf(p) 
for the hardcore potential are plotted in Fig. [h] for 
n = 20. The main features are the following. In the range 
where |whf(p)| is largest, up (p) / uh(p) ~ n 2 ] the abso- 
lute minimum of — uf(p) is constant (given the rescaling 
which keeps vn constant) whereas that of un(p) vanishes 
as n~3. Hence, for large n, the Hartree potential is in- 
significant. 

Following the analysis by FKS, we expect to find a 
periodic sequence of fully occupied and empty stripes of 
Landau orbital guiding centres, the periodicity of which is 
given by the prominent first minimum at p m m of — up(p). 
With p m i n ~ 2.16/Rl, the wavelength is approximately 
A m ; n ~ 2.9Rl, about ten percent larger than that found 
by FKS for the potential Vag- The origin of this differ- 
ence lies in the fact that our contact interaction is less 
hostile to density variations on large lengthscales than a 
Coulomb interaction. 

Higher harmonics do not destabilise this CDW, be- 
cause the oscillations of uf(p) are incommensurate with 
p m in, and because the principal minimum is very pro- 
nounced. 



C. Ground state in a system of infinite size 



As a restricted illustration of this difficulty, we con- 
sider fluctuations starting from a stripe state. Referring 
to Fig. [ll], there are non-vanishing matrix elements con- 
necting the stripe state to a state slightly higher in energy 
with two electrons at the edges moved by an infitesimal 
distance (d < Af) into the empty regions. For a system 
on a finite cylinder, these matrix elements are absent 
due to the discrete locations for the centres of Landau 
orbitals. For the infinite system, however, this problem 
persists, and we address it as follows. The self-energy 
fl5.3| ) will not be changed noticeably by changes in or- 
bital occupation at the edges of the stripes because their 
width is only a fraction n" 1 of the total stripe width, 
and Sfe varies slowly across the edge. After linearising 
Sfe around the edge, the excitations of Fig. [ll] acquire 
a linear dispersion relation. The problem of one isolated 
stripe then corresponds to a Luttinger liquid [ ^2[ p3| , 
where one edge of the stripe provides the left-moving 
fermions and the other, the right moving ones. As a re- 
sult, sharp edges in orbital occupation are smoothed out 
[ E3j, but only over a fraction of the stripe width that 
vanishes for large n. It seems likely that this feature of 
the behaviour should persist for a system of many stripes 
and hence that quantum fluctuations around the stripe 
state are unimportant for a contact potential. 



A detailed study of Hartree Fock ground states for a 
system of electrons interacting via the screened Coulomb 
potential Vag has been carried out by FKS. Whilst it is 
likely that their conclusions apply to a system with the 
hardcore potential, we are unable to make exact state- 
ments either about the lowest energy Hartree Fock state 
or about the importance of quantum fluctuations around 
a Hartree Fock state at zero temperature. 

The Hartree Fock proble m is simplified in princi- 
ple by the result of section VA: we need to find the 



Fourier transform of the electron density, n(q), that min- 
imises J GpqVrff (q)|n(q)| 2 , subject to the constraint that 
n(q) arises from occupation only of states within the 
n th Landau level. In fact, this constraint is very hard 
to incorporate into a rigorous approach. Despite these 
obstacles, FKS, studying Vagi have proposed very plau- 
sible candidates for the ground state on the basis of per- 
suasive physical arguments, supported by extensive nu- 
merical calculations. They find a unidirectional CDW 
ground state for v near 1/2, and a triangular CDW for 
v away from 1/2. Our results at finite temperature, de- 
scribed in section |y|, mirror the conclusions of FKS at 
zero temperature. 

In addition to the difficulty of showing rigorously 
whether the states discovered by FKS lie at the global 
minimum of the Hartree Fock energy, there is a second 
problem in an exact treatment. It is to show whether 
quantum fluctuations around the Hartree Fock state are 
important. 



VI. COMPARISON WITH BEHAVIOUR FOR 
FINITE RANGE POTENTIALS 

A. Diagrammatic analysis 

If the two-particle interaction has a finite range, the 
diagrammatic analysis described in appendix |A| does not 
simplify in the same way as for the hardcore potential. 
To demonstrate the essential differences, we consider, as 



an example of a finite range potential, V2 (2/7). For sim- 
plicity, we include a 'neutralizing background' so that 
V 2 (0) = 0. 

Consider first the exchange diagram with k interac- 
tion lines (A in Fig. ||). Its contribution to the self en- 
ergy scales with n as E ; M(0,0] ~ [v 2 R/Rl] , where 
Rl ~ n-2. We show below that diagrams with crossed 
interaction lines (B in Fig. |^) or k — 1 closed fermion 
loops (C in Fig. |2|), make a contribution of the same 
order to the self-energy. Therefore, these diagrams can- 
not be discarded for large n. We emphasize that, having 
projected the Hamiltonian onto a single Landau level, 
the interaction strength sets the only energy scale. As 
previously, we choose this scale to be our unit of energy, 
setting «2-R/Rl = L 

We now show for the case k — 2 that all three dia- 
grams in Fig. |^ have the same n-dependence for large 
n. We denote the integrals over intermediate states by 
Ia,Ib and Ic and obtain 



G 



Ia = 



qdqw(q) 



Ib = J qdq [w{q)} 

Ic = J J qidqiq 2 dq 2 Jo(<?i92)w(<7iMg 2 ) 

where w(q) = V(q)a n (q 2 /2). It can easily be checked 
that the main contribution to these integrals arises from 
the region where q is of order 1. In this region, w(q) ~ 
n~i. Therefore, I a and Ib are of the same order. The 
Bessel function J (9192) m Ic oscillates slowly in this 
region so that Ic, also, makes a contribution of the same 
order in n. 

The analysis can be extended to all k. The crucial 
difference between between finite range and the hardcore 
interactions is the following. While the contribution to 
I A , Ib and Ic arising from the region q ~ 1 does not differ 
in both cases, the dominant contribution for the hardcore 
interaction arises from 9 > 1, and here the dependence 
on n of Ia,Ib and Ic differs. By contrast, for finite 



range potentials, V(q) is falls off for q 3> R 



1, and 



the large q region never contributes. Since this feature 
is common to all finite range potentials, we expect the 
above argument to be applicable generally. In particular, 
for Vag ; assumi ng that divergences arising in integrals 
such as in Eq. (|Aq) are removed by screening for q < 
n _ 2, we reach the same conclusions up to logarithmic 
corrections. Hence, the Hartree-Fock approximation for 
finite range potentials is not exact even as n — > 00, and 
omits diagrams of the same order in n as those retained. 

Nevertheless, as the range of the potential is increased 
from zero, corrections to the Hartree-Fock approximation 
are governed by the small parameter R/l. For the second 
order diagrams, this can be seen as follows. Consider first 
the integrals I a and Ic- The integral Ic is effectively 
cut off by the oscillations of the Bessel function when 
9i<?2 3> 1, whereas the integrand of I a is not cut off until 
q ~ R , so that Ic/Ia is parametrically small in R/l. 
Similarly, the contribution of the region R^ 1 > q > 1 to 
Ib compared with its contribution to I a is parametrically 
small in R/l, and it is from this region that the dominant 
contribution to Ia arises. 

The condition R/l <§; 1 is fulfilled for the screened 
Coulomb interaction, Vag ; f° r weak fields: taking the 
effective range of Vag to be as, as/Z^l/10 in GaAs 
for a magnetic field B~ 70mT. 

One can still ask whether there in fact exist any com- 
petitors to a CDW groundstate even for R/l > 1. We 
consider one such example: the Tao-Thouless state at 
filling fraction v — 1/t, in which one in t Landau orbitals 
is occupied in a regular fashion. It is of interest because 
it is qualitatively different from the CDW states, having 
a uniform charge density, and it is sufficiently simple to 
allow detailed calculation. We find the cohesive energy 
of the Tao-Thouless state scales as im~ 5 . For compari- 
son, the cohesive energies of the CDW groundstates with 



the interactions V2 and Vag scale as to"'. The nature 
of the true groundstate in high Landau levels with in- 
teractions with range comparable to or larger than the 
magnetic length I, appears to be open and may depend 
on the details of the interaction. 

Tao and Thouless calculated the cohesive energy of 
their state in the random-phase approximation. The 
starting point of the self-consistent calculation is as fol- 
lows: it is assumed that the energy for an occupied Lan- 
dau orbital is rj p and for an empty orbital is This 
determines a gap, A = -q^ — r\ v . The RPA series is then 
summed to obtain expressions for r\ v and r\h in terms of 
A. These equations are solved sclf-consistcntly. General- 
izing the work of Tao and Thouless [ [ll| to higher Landau 
levels and to a wavevector-dependent dielectric constant 
yields the following equation 



Vh - V P 



w(q) 2 qdq 



(AH 2 + 2tAw(q))i + At + 2w(q) 



r°° 2 



w(q) 2 qdq 
At + 2w(q) ' 



(6.1) 



We wish to know the scaling of A, and hence of the 
cohesive energy, with n. By considering the behaviour 
of the integrand as n — > 00, assuming A = Aqu" (with 
Ao independent of n), we find a simply by comparing 
powers of n on both sides of ( |6.1| ) . We find that a = —1/2 
for any fixed t. 

Thus, the cohesive energy of the CDW in the Hartree- 
Fock approximation has the same order as the cohesive 
energy of the uniform density Tao-Thouless state. 



B. Influence of the potential range on the structure 
of the charge density waves 

Even within the Hartree Fock approximation, there are 
significant differences in the structure of the charge den- 
sity waves formed with hardcore or finite range interac- 
tions. Near half-filling, FKS find a state of alternating 
incompressible (fully occupied) and empty stripes of Lan- 
dau orbitals for small and moderate n. Asymptotically 
for large n, however, they find that the stripes break up 
into narrower stripes. When the guiding centre density 
of these narrower stripes is coarse-grained, a sinusoidal 
oscillation in density at wavevector p m i n results. The ori- 
gin of the break-up of the stripes lies in the fact that, for 
n — > 00, the energy cost for higher harmonics in guiding 
centre density outweighs the gain due to the fundamen- 
tal one. This happens because it is the Hartree potential 
which sets the energy scale for higher harmonics, while 
the scale for the energy of the fundamental harmonic is 
set by the Fock potential. At large n, the Hartree poten- 
tial is parametrically larger than the Fock potential. The 
higher harmonics are avoided by the break up into nar- 
rower stripes, which establishes a quasisinusoidal density 
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distribution. For the hardcore potential, this does not 
occur; rather, due to the vanishing of uh(p), the energy 
cost for higher harmonics does not grow relative to the 
fundamental one as n — > oo. 

In order to see whether the break-up found by FKS 
is a general feature of potentials with non-zero range, 
we investigate the t wo simple finite range potentials in- 
troduced in Eq. fl2.7[) . In both cases, we find that 
^mil) = Vi(q)a n (q 2 /2) remains of order Vi for q of order 
as n — > oo. For «j F , this is not the case. The inte- 
gral fif P (p) = J qdqVi(q)a n (q 2 /2)J Q (pq) is effectively cut 
off at q = R , and we find for large n 



1 



%(p) RR L 



(6.2) 



Since the fall off of V(q), and hence of the integral for 
uf(p), above a wavevector q ~ is a feature of all 
finite range potentials, we expect the stripes to break up 
when Rl ^ R- If the effective range of the Coulomb 
potential is taken to be the Bohr radius, this condition 
corresponds to r.T 1 = kpas/V^ ^> 1, which is where FKS 
indeed found the destruction of the wide stripes to occur. 

We note that, for finite range potentials, there also 
exist other, quite different routes by which stripes may 
become unstable. For example, for the potential V±, 
uhf(p) has minima of comparable depth at p ~ n _ 2 and 
at p ~ R , Then, it may be more advantageous to build 
a correlated state where correlations on the scale R are 
as important as or more important than correlations on 
scales of Rl- 

Finally, it should be remarked that the transition to 
a charge-density wave will become increasingly hard to 
observe for n — > oo, since the transition temperature, 
being proportional to uf, decreases with increasing n es- 
sentially as n~ 5 . 



VII. SUMMARY 

In this paper we have considered electrons in a high 
Landau level interacting via a hardcore potential and 
found that this is a rare example of a many-body prob- 
lem in more than one dimension which is, at least in part, 
exactly solvable. 

The focus of our work is the extent to which it is pos- 
sible to obtain mathematically controlled conclusions in 
the high Landau level limit. This turns out to be very 
dependent on the form of the two-particle interaction po- 
tential. For the hardcore potential, we have shown that 
Hartree Fock theory is exact and found a transition to 
both unidirectional and triangular charge-density waves 
at finite temperatures. For interactions of finite range R, 
this is not the case. Nonetheless, for short-ranged inter- 
actions, R/l is a small parameter governing corrections 
to the Hartree-Fock approximation. For R/l > 1, we ar- 
gue that there exist uniform density states which may be 



competitors with the charge density waves for the ground 
state. 

In the context of the lowest Landau level, work with 
the hardcore potential has yielded some significant con- 
ceptual simplifications while retaining the essential phys- 
ical features of the problem [ || . Whether this is also the 
case in high Landau levels will ultimately have to be set- 
tled by experiment. 
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APPENDIX A: BOUNDS ON THE 
CONTRIBUTION FROM NON-HARTREE FOCK 
DIAGRAMS 

We compare the n-dependence of the three types of 
diagrams illustrated in Fig. || for each order k of the 
interaction. For the first class of diagrams (A in Fig. ||), 
it follows from Eq. (|3.5|) that their n-dependence is 



The at are independent of q\ y , q 2y , so that we can use 
the Cauchy-Schwarz inequality to obtain: 



(vn) k 



(Al) 



In the second class of diagrams (B in Fig. §), those 
with crossed interaction lines, we use an analysis follow- 
ing work by Benedict and Chalker [ [24j] on the scattering 
of electrons in high Landau levels with disorder. 

Consider the diagram 

with two crossed lines (Fig. |l^) which corresponds to 
the integral / / dqi y dq 2y M(q 2y ,q ly )M(-qi y ,q 2y ). From 
this, we can guess the general structure for the integral 
over intermediate states I x for r interaction lines. It is 

4 = jY[dqiyM(Qi,q iy ), (A2) 

i— 1 

where Qi = X)j=i ^ij^jy with Zij — if the interaction 
lines i and j do not cross and + 1(— 1) if they cross and i 
starts to the left (right) of j. 

We obtain an upper bound on \I X \ by first picking a 
pair t, u such that Z tu = 1 and then relabeling indices 
such that t = 1 and u = 2. Then 

|4| < I f f[dq jy I 12 M(Qj,q jy )\ 

Il2 = y II d 1*y M (°1 + l2y,qiy) M (a 2 - qiy, q 2 y) 



4 2 | < 



dqi y dq2y\M(qiy,q 2 y)Y 



(A3) 



|4| < Jj\_dq iy \M(Qi,q iy )\ J J dq ly dq 2y \M(q ly ,q 2y )\ 2 (A4) 

i=3 

|4| < [f[dq iy \M(0,q iy )\ J J dq ly dq 2y \M(q ly ,q 2y )\ 2 (A5) 

The last line follows since |M(0,/)| < \M(m,l)\ for 
the hard-core potential. The integrals over qi y (i > 
3) are now decoupled and give the same contribution 
as ( |Al[ ). The remaining n-dependence is fully contained 
in |4 2 | < (2tt)- 2 J qdq(V(q)a n {q 2 /2)) 2 . For our poten- 
tial, (V(q)a n (q 2 /2)) is monotonically increasing up to 
the turning point at q = y8n where it is proportional to 
(vn)n _ 3 . Replacing one of the factors (V(q)a n (q 2 /2)) in 
42 by t his u pper bound, the remaining integral is given 
by Eq. (3.5), so that 42) and hence I x , vanishes at least 
as n~3 compared to the exchange diagram with the same 
number of interaction lines. Therefore, all diagrams con- 
taining crossed interaction lines and no fermion loops 
make a vanishing contribution asruoo. 

Next, consider the class of diagrams containing closed 
fermion loops (C in Fig. ||). First examine those which 
contain Hartree loops (Fig. ||) only. A Hartree loop 
gives rise to an integration of the form ^ m M(m, 0) 
= (2tt)- 2 / / dq x dmV{q x )a n {ql/2)e- i 'i- m = 0, since 
Vhc(0) = 0. Similarly, all other diagrams containing 
Hartree loops vanish for our potential, a fact unrelated 
to the large-n limit. Note that this is no longer the case 
in a state with broken translational symmetry. 

Finally, we turn to those diagrams which contain at 
least one non-Hartree fermion loop such as the one in 
Fig. |l3|. Each of those loops give rise to an integral 
over the loop pseudomomentum s, which leads to a 8- 
function 8(qi x + q 2x — qz x ). Similarly, pseudomomentum 
conservation gives rise to a corresponding condition on 
the qiy. The integral over intermediate states, 4 of di- 
agrams containing k interaction lines and L non-Hartree 
loops is then bounded above by 



k-L 



141 < 



:I1 / ^XIIMGOVW)! 



(A6) 



1=1 



where Qi = X^=i ^jQj' an d Yij — 0, ±1 is fixed by the 
constraints imposed by the <5-functions. This is an upper 
bound since we have left out the oscillatory term in the 
expression for the vertex. 

We can again replace the \ot n (Qi)V(Qi)\ by an upper 
bound and finally obtain ~ (vn) fe (n~3) L ~ n~~ . 
Hence these diagrams also make a vanishing contribution 
asm 00. 



9 



APPENDIX B: CALCULATION OF THE FOCK 
POTENTIAL 

We start by writing the Fourier transform of the real 
space density, n(q), in terms of the density matrix p kt y in 
the Landau orbital basis. 

ft(q) - X)< fc/ I e<q "l*)^*' 
k,k' 

L 

~ 2tt 



dke ik ^p ktk+qy ) exp(^K(|- ) 



= ^P(fe,ft)exp(— ^— )w„(y) 



(Bl) 



This can be inverted to give p(q x , q y ), and hence p k>k i , 
in terms of n(q). 

The Fock potential in terms of the real space density 
matrix p(r, r') is 

y F (r,r') = V(r-r>(r,r'). 



(/|y F |m> = £ | AP (/| e ^-|fe) 



(B2) 



Substituting for /j^/ from ( |Bl|) gives 
F F (r, r') = J2 V(r - r')(r'|fc)(fc + q y \v)p k , k+qy (B3) 

and 

J2 T 
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FIG. 3. Hartree loop. The solid line is a full propagator. 
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FIG. 4. Dyson equation for the self-energy. 
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FIG. 5. The self energy for the high-temperature phase. 



(k 



\m)p ktk+qy V{p) 



d 2 q . 



(2tt) 2 

where the Fock potential is 
t?f f s 1 f d 2 P T ~ 



n(q)(/|e-^ r |m)^ ff (q), 



V(p)e- tp - q a„(i-). 



(B4) 
(B5) 

(B6) 
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FIG. 1. Interaction vertex M(m,l). 



FIG. 6. Ladder diagrams for the two-particle propagator. 
The solid lines denote full propagators. 
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FIG. 7. Bethe-Salpeter equation for I\ Greek letters de- 
note Matsubara frequencies. 



FIG. 2. The three classes of diagrams. 
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FIG. 8. Phases near the triple point. The solid line, Ti, 
separates the uniform and triangular phases, the short dashes, 
T2, are the spinodal line and the long dashes, T12, separate 
the triangular and unidirectional phase. 
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FIG. 10. uh(p) and — uf(p)- The unit of p is the inverse 
magnetic length; the unit of energy is arbitrary. 
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FIG. 9. Phase diagram as in the previous figure. The dot- 
ted line is an extrapolation of T12. 



FIG. 11. Two near degenerate states \A), \B) with nonva- 
nishing matrix element (A|Vhc|-B). 
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FIG. 12. Diagram with two crossed interaction lines 
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FIG. 13. Closed fermion loop. 
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